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Abstract: We use holography to study J\f = 4 supersymmetric SU{Nc) Yang-Mills the- 
ory in the large- A^c and large-coupling limits coupled to a number Nf <^ Nc of (n -|- 1)- 
dimensional massless supersymmetric hypermultiplets in the fundamental representation 
of SU{Nc), with n = 2,3. We introduce a temperature T, a baryon number chemical po- 
tential n, and a baryon number magnetic field B, and work in a regime with n ^ T, ^Tb. 
We study the collective excitations of these holographic quantum liquids by computing the 
poles in the retarded Green's function of the baryon number charge density operator and 
the associated peaks in the spectral function. We focus on the evolution of the collective 
excitations as we increase the frequency relative to T, i.e. the hydrodynamic/collisionless 
crossover. We find that for all i?, at low frequencies the tallest peak in the spectral function 
is associated with hydrodynamic charge diffusion. At high frequencies the tallest peak is 
associated with a sound mode similar to the zero sound mode in the collisionless regime 
of a Landau Fermi liquid. The sound mode has a gap proportional to and as a result 
for intermediate frequencies and for B sufficiently large compared to T the spectral weight 
is strongly suppressed. We find that the hydrodynamic/collisionless crossover occurs at a 
frequency that is approximately i?-independent. 
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1 Introduction 

Consider a system with a global U (1) symmetry, in a state with a net charge density of that 
C/(l). Such a state is compressible if the charge density varies smoothly as a function of 
the chemical potential. The best-understood examples of compressible states are solids, in 
which translational symmetry is spontaneously broken to a discrete subgroup, superfluids, 
in which the U{1) is spontaneously broken, and Landau Fermi Liquids (LFLs), in which 
neither translational symmetry nor the U{1) are broken. 

The degrees of freedom in a LFL are long-lived fermionic quasi-particle excitations 
about a Fermi surface. Two characteristic features of a LFL are that at temperatures T 
low compared to the Fermi energy, the heat capacity cy oc T and the electrical resistivity 
p ocT^. The spectrum of a LFL includes not only particle- hole excitations about the Fermi 
surface but also collective charge density excitations, which produce poles in the retarded 
Green's function of the charge density operator. At sufficiently low T, in the so-called 
"collisionless" regime in which quantum effects are more important than thermal effects, 
the spectrum of collective excitations includes "zero sound," a fluctuation of the shape of 
the Fermi surface that has the dispersion relation of a sound wave [1-3]- At sufficiently high 
T, in the hydrodynamic regime in which thermal effects are more important than quantum 
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effects, the spectrum of collective excitations includes the expected hydrodynamic sound 
and charge diffusion modes [3]. The transition between these two regimes as T increases 
is called the collisionless/hydrodynamic crossover. 

Some real materials are compressible, break neither translational symmetry nor the 
C/(l), and exhibit Fermi surfaces, but the low-energy degrees of freedom are not those of a 
LFL. These materials are generically called "non- Fermi liquids" (NFLs). A subset of these 
are the so-called "strange metals," including the normal (non-superconducting) phase of 
high-Tc superconductors. The most salient characteristic feature of a strange metal is a 
linear resistivity, p (x T. Experimental evidence indicates that, in many cases, the strange 
metal phase is related to quantum criticality, i.e. a continuous phase transition at zero 
temperature [4]. In particular, the degrees of freedom at such a quantum critical point, 
when heated up to non-zero T, may give rise to a quantum critical region in which p (xT. 
Demonstrating this from first principles is difficult because generically quantum critical 
degrees of freedom are strongly interacting. The existence of strange metals and other 
NFLs raises a crucial question: can we classify compressible ground states? 

The Anti-de Sitter /Conformal Field Theory (AdS/CFT) correspondence [5-7], and 
more generally holography, may help answer this question. In AdS/CFT the conserved 
current of a global ^7(1) symmetry of the CFT is dual to a U{1) gauge field in AdS, 
and compressible states can be described by asymptotically AdS spacetimes with non- 
zero electric fiux at the boundary. Holography provides many examples of compressible 
states that are not solids, superfiuids, or LFLs, the canonical example being the planar 
Reissner-Nordstrom-AdS solution of Einstein-Maxwell theory [8]. Holography has been 
used to describe both quantum criticality [9-16] and compressible states with p oc T [17- 
19] and, more generally, has the potential to reveal some guiding principles for classifying 
compressible states [20, 21]. 

With these motivations in mind, we will study two holographic systems describing 
compressible states that are not solids, superfiuids, LFLs, or NFLs, and which exhibit 
quantum phase transitions for sufficiently large magnetic field B. Our main goal is to study 
the effect of non-zero B on the collisionless/hydrodynamic crossover in these systems. 

We will study the field theories arising from the (n + l)-dimensional intersection of Nc 
D3-branes with A'^^ Dp-branes in type IIB string theory, with p = 2n + l = 5,7so that 
n = 2, 3 [22-25]. We will call these the D3/Dp systems or D3/Dp theories. These theories 
are (3 + l)-dimensional = 4 supersymmetric SU{Nc) Yang-Mills (SYM) theory coupled 
to a number Nf of (n + l)-dimensional hypermultiplets in the fundamental representation 
of SU{Nc), i.e. fiavor fields. We consider only massless fiavor fields. We will work in the 
't Hooft large- limit, with large 't Hooft coupling, and in the probe limit Nf <C Nc- 

Each D3/D|3 theory enjoys a global U{l)f, baryon number symmetry, with conserved 
current J'^, where i' = t,x,y and for n = 3 also u = z. To study compressible states of these 
theories, we will introduce a baryon number charge density (J*), or equivalently a baryon 
number chemical potential p. Recall that a hypermultiplet contains fermions and scalars, 
and that both are charged under U{l)b. States with non-zero (J*) are dual holographically 
to probe Dp-branes in AdS with non-zero worldvolume electric flux [22, 23, 26]. 

For temperatures T <C these compressible states have been studied extensively, using 
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holography, in refs. [27-36]. Neither translational symmetry nor U{l)i, are broken in these 
states, and to date no evidence of a Fermi surface in these systems has been found. In short, 
the evidence accumulated so far suggests that these states are not solids, superfluids, LFLs, 
or NFLs. Indeed, these states exhibit various unusual properties, including an extensive 
T = ground state degeneracy [28] . We will call these low-temperature compressible states 
"holographic quantum liquids." 

The spectrum of excitations in these holographic quantum liquids at T = includes 
sound modes, i.e. two modes producing poles in J*'s retarded Green's function, G^}^{uj,k), 
that, when expressed as complex frequencies a;, have real parts proportional to momentum 
k and imaginary parts oc k"^ [28]. Due to the similarity with the zero sound mode in a LFL, 
we will call these modes "holographic zero sound." 

As T increases, these holographic quantum liquids exhibit a collisionless/hydrodynamic 
crossover analogous to that in a LFL. For the D3/D7 system the crossover was studied in 
ref. [33], and a similar crossover occurs in the D3/D5 system, as we will show. The crossover 
is simplest to define from the behavior of poles of G^(a;,A;). We can identify poles as 
holographic zero sound or charge diffusion by computing their dispersion relations. When 
k S> T^/ 11, we can identify the two poles of G''^{uj, k) closest to the origin of the complex uj 
plane as holographic zero sound. As T increases, these poles move down into the complex 
uj plane and towards the imaginary axis. When k ~ T"^ j \i, the two poles collide on the 
imaginary axis. As T continues to increase the poles split, one pole moving up and one 
moving down the imaginary axis (see figure 5). When k <C TVa*, we can identify the pole 
moving up the imaginary axis as that of charge diffusion. We define the crossover as the 
point where the poles collide on the imaginary axis. The crossover is also apparent in J*'s 
spectral function, which as a function of a; exhibits a single peak due to holographic zero 
sound when k 3> and a single peak due to charge diffusion when k <C T"^ I y^. 

Now let us introduce a non-zero f/(l)t, magnetic field B, which appears in the bulk as 
magnetic flux on the Dp-brane worldvolume [37]. When T = and (J*) is non-zero, in 
each of the D3/Dp systems a sufficiently large B triggers a quantum phase transition in 
which a global symmetry is spontaneously broken [9-12]. For n = 3 the transition is at 
least second order [10], while for n = 2 the transition is of BKT-type [11, 12]. When T is 
non-zero, each quantum critical point gives rise to a quantum critical region. 

We want to determine how B affects the collisionless/hydrodynamic crossover in these 
holographic quantum liquids. Five scales appear in our problem: /x, T, B, which character- 
ize the state of the system, and u, k, which characterize perturbations of the state. We will 
work in the quantum liquid regime, where fi ^ T, \fB. Furthermore, we will always work 
with values of B below any known phase transition, where we already find a rich story. We 
choose the momentum to point in the x direction. For n = 3 we choose the magnetic field 
to point in the z direction, orthogonal to the momentum. We work with w. A; <C //. Our 
results are qualitatively similar for both of the D3/Dp theories. 

In our systems, a non-zero B has two major effects on the spectrum of collective 
excitations at low temperatures. First, a non-zero B breaks parity, and so operators can 
mix that could not mix when i? = 0. In particular, we find that when B is non-zero the 
poles in G^l^{uj,k) mix with those of J^'s retarded Green's function, G^^{uj,k). Due to 
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the mixing of poles, for small T the holographic zero sound poles are no longer the only 
poles near the origin of the complex u plane: a purely imaginary pole also appears. As a 
result, the spectral function of J*, as a function of uj, can exhibit two large peaks, one from 
holographic zero sound and one from the purely imaginary mode. 

Second, when T = 0, a non-zero B produces a gap in the holographic zero sound 
dispersion relation (x B/fi [38, 39]. Such a gap is consistent with Kohn's theorem [40], 
which states that for non-relativistic particles with pair-wise interactions of any strength, 
in an external magnetic field the center-of-mass motion decouples from the relative motion 
of the particles and behaves as a free particle in a magnetic field, which then guarantees a 
k = resonance at the cyclotron frequency, which is B. 

For sufficiently small B, meaning B <T^, and for fixed k, we find that as T increases 
a collision of poles on the imaginary axis still takes place. Due to the mixing of poles, the 
details of the collision differ from the B = case, nevertheless we can still use the collision 
of poles to define the location of the collisionless/hydrodynamic crossover. For sufficiently 
large B, meaning B >T^, however, no such collision of poles occurs, so at first glance we 
seem to be unable to identify a crossover at a specific value of T. 

After a thorough search, however, we have found a simple alternative definition that 
allows us to identify a precise location for the crossover, for all values of B that we study. 
If we express the poles in G^^{uj, k) in terms of real lo and complex k, then in the complex 
k plane we always find a single pole near the origin, for all values of B that we study. For 
low T, meaning to ^ T'^ / fJ-, and for lo larger than any other scale except //, we can identify 
the pole as that of holographic zero sound. As T increases, the pole becomes that of charge 
diffusion when oj ~ 0.45 T^/^u for the n = 3 case and when oj ~ 0.30 T^/^u for n = 2. Most 
importantly, these values of T are independent of B, and indeed are exactly the same values 
as those at 5 = 0. We thus define the location of the collisionless/hydrodynamic crossover 
as the point where the pole in the complex k plane becomes that of charge diffusion. 

What happens to the pole in the complex k plane between lo ^ T'^/fi and co » T^/fJ- 
depends on the value of B, and is simplest to understand by studying the spectral function 
of J*. For any B, we find that the spectral function, as a function of /c at a fixed uj, exhibits 
only a single peak. For any B, when ui < r^//i the peak is due to charge diffusion. For 
small B, as we increase to, the peak changes to that of holographic zero sound at a; 
At large B, however, when w w T"^ / 1^ the peak is suppressed by orders of magnitude, 
and only when oj k, B / ^ does the peak grow again, becoming that of holographic zero 
sound. The suppression is simple to understand: once B >T^ the holographic zero sound 
dispersion relation acquires a gap oc B/n, so if we perturb the system with a frequency 
in the collisionless regime, u > T'^ / fJ^, but with oj < B / fi, then the frequency is not high 
enough to overcome the gap and excite the holographic zero sound. We thus find an 
intermediate region T^//i ^ <^ < B/^, with relatively low spectral weight. We summarize 
these results schematically in fig. 1 (see also fig. 22). 

We have several reasons to believe that results qualitatively similar to ours will appear 
for any compressible state with a holographic description in terms of probe D-branes with 
worldvolume electric fiux. The only ingredient in our calculation is a probe D-brane in AdS 
space. We see qualitatively similar results for both of our D3/Dp systems. A gap appears 
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Figure 1. A schematic summary of our main results. We depict the plane of B/T^ versus ujfi/T'^, 
with colors indicating the nature of the peak in J*'s spectral function, at fixed w, as a function of k. 
In region (I), shaded red, the peak is due to charge diffusion. In region (III), shaded blue, the peak 
is due to holographic zero sound. In region (II), shaded grey, the peak is not associated with either 
charge diffusion or holographic zero sound, and in fact in region (II) the peak is suppressed by 
several orders of magnitude compared to the peaks in regions (I) and (III). The suppression is due 
to the gap oc B/fi in the holographic zero sound's dispersion relation, which defines the diagonal 
line between regions (II) and (III). We define the coUisionless/hydrodynamic crossover as the point 
where, as we increase T, the pole in J*'s retarded Green's function becomes that of charge diffusion. 
The crossover so defined is independent of B, and appears at the vertical line at the right-hand 
boundary of region (I). 



in the holographic zero sound dispersion relation in several different holographic probe 
D-brane systems [38, 39, 41]. Indeed, Kohn's theorem is very general, and suggests that at 
high frequencies and high magnetic fields, compared to T, we should always find a region of 
suppressed spectral weight, like region (II) in fig. 1. In short, our fig. 1 may be characteristic 
of holographic quantum liquids described by probe D-branes with worldvolume electric flux, 
and hence may help in classifying them. 

This paper is organized as follows. In section 2 we review the thermodynamics of our 
holographic quantum liquids and in section 3 we review how to compute poles in retarded 
Green's functions, and spectral functions, using holography. In section 4 we discuss the 
collisionless/hydrodynamic crossover in our systems when B = 0, which for the D3/D7 
system is a review but for the D3/D5 system includes novel results. In section 5 we discuss 
the crossover in our systems at non-zero B. We conclude in section 6 with some suggestions 
for future research in these (and similar) systems. 



-5- 



2 Review of Holographic Quantum Liquids 



In type IIB string theory, we will study the (n+l)-dimensional intersection of Nc D3-branes 
with Nf Dp-branes, with p = 2n + 1 = 5, 7, which we call the D3/Dp systems [22-25]. We 
summarize both systems in the following array: 
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In these intersections the low-energy theory on the D3-brane worldvolume is (3 + 1)- 
dimensional = 4 SYM theory with gauge group SU (Nc) and Yang- Mills coupling squared 
9ym ~ with string coupling Qs, coupled to A'^^ {n + l)-dimensional massless hyper- 

multiplets in the Nc representation of SU{Nc), i.e. flavor fields, preserving eight Poincare 
supercharges. In the n = 2 case the flavor fields propagate along a codimension-one de- 
fect, which we take to be the plane z = 0. We will take the 't Hooft limit Nc ^ oo with 
9ym ~^ 0, with the 't Hooft coupling A = gyj^^Nc fixed, followed by the large-coupling limit 
A ^ 1. We will always work in the probe limit: we keep Nf fixed as Nc ^ oo, expand all 
observables in the small parameter Nf/Nc, and retain terms only up to order NfNc. 

The R-symmetry of = 4 SYM is SO{Q), corresponding to rotations in the (X^, . . . , X^) 
directions. The (n + l)-dimensional fiavor fields break the SO{6) to SO{n + l) x S0{5 — n), 
corresponding to rotations in the {X^, . . . ,X'^~^^) and (X"^^, . . . ,X^) directions, respec- 
tively. In the n = 2 case the S0{3) x S0{3) is the R-symmetry of the remaining (2 + 1)- 
dimensional A/" = 4 supersymmetry. In the n = 3 case an SU{2) x U{1) subgroup of 
50(4) X U{1) is the R-symmetry of the remaining N = 2 supersymmetry. 

We want to study thermal equilibrium states of these theories at temperature T. A 
thermal state of large- A^o strongly-coupled AA = 4 SYM is dual to type IIB supergravity 
in the near- horizon geometry of non-extremal D3-branes, j4(iS5-Schwarzschild times an S^, 
whose metric takes the form 

ds^ = gtt{r)dt^ + g.^.^{r) {dx^ + dy"^ + dz^) + grr{r)dr^ + R^ds% 

= -^f(^)df + ^ {dx' + V + dz') + R'ds%, (2.2) 

where 

m = i-"fl. (2.3) 

In eqs. (2.2) and (2.3), r is the AdS^ radial coordinate, with the AdS^ boundary at r — )• c« 
and the planar Schwarzschild horizon at r = rn, R is the AdS^ radius of curvature, which 
is fixed to be i?^ = AirgsNca''^ , with a' the string length squared, and ds^s is the metric 
of a unit-radius . The full type IIB supergravity solution includes Nc units of Ramond- 
Ramond (RR) five- form flux on the . The Hawking temperature T of this planar black 
hole is given by = -kTE? . We also identify T as the temperature of the dual field theory 
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state [42]. Starting now, we use units in which R = 1. In these units we can convert 
between supergravity and field theory quantities using a' = 1/ \/A. 

The massless probe flavors are dual to probe Dp-branes extended along AdSn+2 x S" 
inside the ^(i55-Schwarzschild times S^, where the AdSn+2 is spanned by {t,x,y,r) and 
for n = 3 also z, and the 5" G [22, 23]. The 5^ has an 50(6) isometry, dual to the 
50(6) R-symmetry of M = i SYM. The Dp-brane breaks that to SO{n + 1) x 50(5 - n), 
dual to the corresponding global symmetries in the field theory. The SO{n + 1) factor is 
the isometry of the 5". 

The non-Abelian Dp-brane action is not known in full generality, however in what 
follows we will only need the Abelian Dp-brane action, which consists of two types of terms, 
a Dirac-Born-Infeld (DBI) term and Wess-Zumino (WZ) terms describing the coupling to 
background RR fields. In what follows the WZ terms will play no role, so we will omit 
their explicit forms. The action for our Dp-branes, 5dp, is thus the DBI term alone. 



where Top = (27r)~PgJ^a'~(P"'"^)/^ is the Dp-brane tension, ^'^ with a, 6 = 0, 1, . . . ,p are the 
Dp-brane worldvolume coordinates, g^^ is the induced metric of the Dp-brane, and Fab is 
the field strength of the U{1) gauge field Aa on the Dp-brane.^ 

The (9-p) scalars propagating on the worldvolume of the Dp-brane, which describe the 
embedding of the Dp-brane in the background, enter 5dp via the induced metric . In 
our cases, of the (9-p) scalars, the 5 — n describing the position of the Dp-brane on the 5^ 
form a vector of SO {5 — n), one of which is holographically dual to the supersymmetric 
mass operator Om- We will always work with trivial worldvolume scalars,^ and as a gauge 
choice we will identify the worldvolume coordinates with the background coordinates along 
the AdSn+2 X S"' spanned by the Dp-brane. In that case, g^^ coincides with the background 
metric in the AdSn+2 x subspace. 

The conserved U{l)h current J'^ is dual to A^,. We want to study states with non- 
zero (J*) and B. We will demand that these states preserve various symmetries: time 
translations, n-dimensional translations, 50(n-|-l) x 50(5— n), and the reflection symmetry 
z —7- —z. To describe non-zero (J*) and preserve these symmetries, we introduce At{r). 
Employing Ar = gauge, we thus have a worldvolume electric flux Frt{r) = A'i[r)^ where 
prime denotes dr- To describe non-zero B and preserve these symmetries, we introduce 
the magnetic flux F^y = B, which must be constant to satisfy the Bianchi identity. With 
these fluxes and with trivial worldvolume scalars, the Dp-brane action reduces to 



^Notice that with respect to the usual conventions (such as those of ref. [43]) we have absorbed a factor 
of (27ra') into the field strength. 

^The one and only place where we discuss solutions with non-trivial worldvolume scalars is at the end 
of this subsection, when we briefly review the solutions of refs. [9-12]. 
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where Np = N fT-QpVs-n- ■, with Vgn the volume of a unit-radius S"". Exphcitly, 



Mp ^ NfT^pVs^ = {WJ J _ ; (2-6) 



(2^) 



NfNc p = 7, 



where in the second equahty we have converted to field theory quantities. In eq. (2.5), 
V^(„,i) is the (infinite) volume of M*^"''"'^^. Starting now, we will divide both sides of eq. (2.5) 
by V^(n,i) ) and in an abuse of notation we will use S^p to denote the resulting action 
density, 5'Dp/^(n,i) — ^ S^p, which we will henceforth refer to as the Dp-brane action. 

Since 5'dp depends only on the derivative of At{r), we obtain a constant of integration, 
-j^, which as shown in refs. [26, 27] determines the charge density in the field theory: 

(J*) = (2W)^ = AA,,(r^)/^7^i:T^^2fi^. (2.7) 



Solving for A'^, we find 



where 



A[ir) = ^ '^V^^ (2.^ 

'g^^' (gL + B') + <P 



We will use "charge density" to refer to either (J*) or d. 

We will now review various properties that characterize the compressible states de- 
scribed by the solution in eq. (2.8), as determined from holographic calculations. We begin 
with low but non-zero temperatures, T <C d^^", and with B = 0. In this regime, in 
these compressible states neither n-dimensional translations nor U{l)i, are broken. In the 
canonical ensemble, the chemical potential is [28] 

^(T, d) = ^ ~ fl ^ - i^T) [1 + (T^Vd^)] . (2.10) 

1 (1/2) 

From eq. (2.10), we see that when T = the density d^/" oc //, so the density varies 
smoothly as a function of fj,, i.e. these states are indeed compressible, and a large density 
implies a large chemical potential. For convenience, we will denote the value of the chemical 
potential at T = as 

— fm — ^ • (^-^^^ 

At low temperatures, the leading density-dependent contribution to the heat capacity is 
cy cx T^"/(i [28], which is markedly different from that of a LFL, where cy cx T, or a gas 
of free bosons, where cy oc T", while the leading contribution to the U{l)b resistivity is 
/9 oc [44], the same as in a LFL. 

Next let us consider T = 0, still with B = 0. The system actually has a non-zero 
T = entropy density s oc (J*)/\/A [28], indicating some degeneracy of states, in stark 
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contrast to solids, superfluids, and LFLs, which are all unique ground states. Degeneracy 
suggests instability, since generically any perturbation will break the degeneracy and drive 
the system to a new, presumably unique, ground state. These states are known to be 
thermodynamically stable, in the sense that they are at least local minima of the free 
energy [45]. As shown for the n = 3 case in ref. [32], these states are also stable against 
dynamical (non-zero oj and k) fluctuations. To date, no Fermi surfaces have been detected 
in these states [10, 32, 38]. Remarkably, these states are in fact merely isolated points in 
an entire moduli space of compressible states parameterized by the expectation values of 
certain scalar operators, as discussed in detail in refs. [35, 36]. The existence of the moduli 
space may be related to the non-zero T = entropy density [36] . 

All of the above facts together indicate that these states are not solids, superfluids, 
LFLs, or NFLs. Indeed, as proposed in ref. [28], these states may be examples of some new 
kind(s) of compressible matter. Some evidence even suggests that the low-energy theory 
describing excitations about these states may be a (0+l)-dimensional CFT [11, 31, 32, 38]. 

Now let us turn to the thermodynamics of these states with non-zero B. For the n = 3 
case, a non-zero B explicitly breaks parity and charge conjugation symmetries, and breaks 
rotational symmetry from 50(3) down to SO{2). Moreover, with T = and non-zero 
(J*), a sufficiently large B triggers a quantum phase transition [9, 10]. When B is below 
a critical value, B < Be, in the ground state (J*) is non-zero but (Om) is zero. When 
B > Be, the ground state changes such that (J*) and (Om) are both nonzero. In the bulk, 
the worldvolume scalar dual to Om becomes nontrivial when B > Be- The non-zero {Om) 
signals spontaneous breaking of the SU{2) x U{1) symmetry down to SU{2). The system 
thus exhibits a quantum phase transition. Notice that U{l)h is unbroken for all B. A 
numerical analysis indicates that the critical value of B is Be ^ 2.138 cP^^ and that the 
transition is not first order, i.e. is at least second order [10]. At non-zero T, the quantum 
critical point gives rise to a quantum critical region [9, 10]. 

For the n = 2 case, a non-zero B breaks parity and charge conjugation symmetries. 
Here again, with T = and non-zero (J*), a sufficiently large B triggers a quantum 
phase transition [11, 12]: for B < Be (J*) is non-zero and {Om) = 0, while for B > 
Be both (J*) and {O m) are non-zero, and the non-zero {Om) spontaneously breaks the 
50(3) X SO{3) symmetry down to 50(3) x U{1). In the bulk, the worldvolume scalar 
dual to Om becomes nontrivial when B > Be- The critical value of B is known exactly. 
Be = ~ 0.378 d, and the transition is of BKT-type: schematically, as B increases 

past Be, {Om) oc e^V'^/-^'=~^/-^ [11]. At any non-zero T the transition becomes second 
order, and the quantum critical point gives rise to a quantum critical region [11, 12]. 

We hasten to add that we have only reviewed the currently-known phase diagrams of 
these systems. For given values of T, B, and d, other states may exist with lower free 
energy than any of the states mentioned above. Indeed, in the n = 3 case, when B is 
non-zero an instability has been detected at a single point in the grand canonical phase 
diagram, (T, //) = {^y/B,y/B) [46], where the known ground state has non-zero (J*) and 
{Om) = 0. The instability has non-zero k, suggesting the existence of an inhomogeneous 
state with lower free energy than the known ground state [46]. In what follows, we always 
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work in regions of the phase diagram where the known ground state has non-zero (J*) and 
{Om) = 0. In bulk terms, we exclusively consider values of T, B and d for which the only 
known solutions for the Dp-brane world volume scalars are the trivial solutions. 

3 Spectral Functions from Holography 

To probe a many-body system experimentally, we can perturb the system by a source for 
some operator. The simplest perturbations to study have either a fixed frequency or a 
fixed momentum. Within linear response theory the change in the expectation value of an 
operator, such as our U{l)b current J^, is given in terms of the small-amplitude source 

au{uj,k) by 



where G^{uj,k) is the retarded Green's function. In general, G^{u,k) is a complex 
function of the two real variables, cj and k. If we naively complexify either uj or k, then 
often poles appear in G^{u}, k), indicating a large response to an infinitesimal source. The 
position of each pole can be expressed as a function of the complex variable in terms of 
the real variable (as well as T, B, etc.); this complex function is called the dispersion 
relation of the excitation corresponding to the pole. Henceforth we will abuse notation and 
write the dispersion relations simply as to (k) or k{u}). Using the Fourier decomposition, 
we can interpret Reu;(A;) as the propagating frequency and —lmuj{k) as the decay rate of 
the excitation, or alternatively Iiek{uj) as the propagating momentum and ImA;(a;) as the 
attenuation, or decay length, of the excitation. 

The poles described above are not directly observable: in a real experiment, neither 
u nor k can be complex. A common experimental technique for studying excitations of a 
many-body system is to measure the spectral function 



which is a real-valued function of real lo and k, and hence is observable. Physically, 
X^'^{u},k) determines the rate of work done on the system by the small external source 
ay{uj,k) [47]. A pole in G^{uj,k) with sufficiently large residue will produce a large peak 
in x^'^i^jk). We can thus identify excitations of the system from peaks in x^'^i^^k). In 
general, the precise relation between the dispersion relation of a pole in G^{io,k) and 
the location of the corresponding peak in x^'^i^jk) is complicated, but typically a pole 
with dispersion relation lo (k) produces a peak in x^'^ i^) centred near Rea;(/c), with width 
proportional to — Ima;(A;). Notice that given x^'^i^^k) we can reconstruct ReGj^{ijj,k) 
using the Kramers-Kronig relation, provided the large-cj and large-A; asymptotics have been 
suitably regulated [47] . 

In holography the on-shell bulk action is the field theory generating functional [6, 7], 
hence to compute a two-point function such as Gj^{u),k), we need to compute the bulk 
on-shell action and take two functional derivatives, which in turn means we need to solve 
the linearized equation of motion for the bulk fields. We thus introduce fluctuations of all 
the Dp-brane worldvolume fields, including all components of Aa and all (9-p) scalars. The 



5 {Jf'ioj, k)) = G^j^{uj, k) a^{uj, k) 



(3.1) 



X' 



'^%u,k)^i G^j^{u:,k)-G^^{u,k)^ , 



(3.2) 
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fluctuations form representations of the symmetry preserved by the background solution. 
At the hnearized level different representations cannot couple. In our case the background 
solution eq. (2.8) preserves SO{n + 1) x SO{b — n) and the reflection symmetry z — t- —z. 
Under S0{n+1), the components of Aa on the 5"" are in the n+1 while all other fluctuations 
are singlets. Under 5*0(5 — n), the worldvolume scalars describing the position of the Dp- 
brane on the are in the 5 — n and all other fluctuations are singlets. We are interested 
only in J*^, which is a singlet under SO{n + 1) x 50(5 — n), and so is dual to a fluctuation 
of Ai, that is also a singlet, and thus has no dependence on the 5" directions and will not 
couple to fluctuations of Aa on the 5" or to the worldvolume scalars describing the position 
of the Dp-brane on the 5^. In the n = 2 case, under the z — )• —z symmetry the fluctuation 
of the worldvolume scalar describing the position of the D5-brane in the z direction is 
clearly odd, and hence does not couple to the a^, which are even. We thus only need to 
consider the fluctuations of A,^, 



On the right-hand-side of eq. (3.3) the first two terms produce the electric and magnetic 
fluxes Frt{r) and Fxy = B of the background solution. We will consider fluctuations a^^ 
that depend on r, t, and x only. The background solution preserves time translations and 
n-dimensional spatial translations, allowing us to perform a Fourier decomposition. 



To preserve Ar = gauge, we take = 0. 

The background solution preserves 50(2) rotational symmetry in the {x, y) plane. 
Under this 50(2), at is a singlet while and ay together form a doublet. Our choice of 
momentum only in x breaks this 50(2). In the D3/D5 system we make that choice with 
no loss of generality, but in the D3/D7 system we are restricting to a special subset of 
fluctuations with momentum perpendicular to the magnetic field. In the D3/D7 system, 
our choice of momentum preserves the z — )■ —z symmetry, under which is odd and hence 
will not couple to at, ax, or ay, which are even. When B = 0, with our choice of momentum 
at and ax can couple to one another, but parity forbids these from coupling to ay. When 
B is non-zero and parity is broken, all three fluctuations may couple to one another, and 
indeed do, as we will see shortly. 

If different components of in the bulk couple, then the dual operators will mix, and 
in particular the poles in different components of G^{uj,k) will mix. These poles will be 
shared by all the components of G^{oj, k) that mix, but the spectral functions can still be 
different since the residues of the poles may be different for each component. 

The quadratic action of the fluctuations, 5^2, and the resulting equations of motion 
are straightforward to obtain but unilluminating, so we will omit them. We will record the 
equation for Or, which in Or = gauge is a constraint on the other fluctuations, 



A^{r, x") = At{r) 5\ + Bx 5\ + a^{r, x^'). 



(3.3) 




(3.4) 



oj at + u{r)'^k a'^ = 0, 



(3.5) 
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where 



u{ry 



\9tt\gr 



A? 



\9tt\gxx 



9r 

Qxx 



(3.6) 



With our choice of momentum, the gauge-invariant fluctuations are ay itself and also the 
electric field (in Fourier space) [48] 



E{r^ uj,k) = k at{r, io,k) + u ax{r, uj, k), 
which is dual to the operator 

= kJ* + uj,r. 



(3.7) 



(3.8) 



In terms of these gauge-invariant fluctuations, using eq. (3.5) and performing an integration- 
by-parts, we can write 5*^2 as 



dr 



dujdk 



\9tt\gxx 



1 .wi2 g_ 

\gtt 



(3.9) 



\E 



^2 _ y:(^j.'^2j^2 

u{r)gU' {gL + B') 



I \2 I g-rr I 2 
S + I 1 



u{r) k 



+iB 



\gu\gxx 



dr 



' (n-l)/2 -1/2 
gxx grr -^t 



where E* and a* denote complex conjugates. For the D7-brane, in eq. (3.9) we have omitted 
a boundary term produced by the integration-by-parts since this boundary term vanishes 
on-shell. The equations of motion that follow from £'^2 are 



E" + 



5rlog 



(n+l)/2| I -1/2 
yxx I ytt I yrr 



+iB 



u[r 



(^2 _ uirfk"^) u{r) {gl^ + B'^) 
){glx + B^)g 



E' + ^ / 2 



\gtt\ 



u{rYk^) E (3.10a) 



dr log 



\gtt\gxx 

' (n+l)/2| 1-1/2^ 

gxx \gtt\grr 



dr 



An-l)/2 

gxx 



ijrr -^t 



-iB 



^ u{r) (5L + B^) 
u{r) [gl, + B^) g 



ay + 



<r) {gl, + 52) 

grr / 2 t \2i 2 



U 



{rfk") 



igtti 



{u^ - u{rYk^) ay 



1/2' 
rr 



\gtt\gxx 



dr 



' (n-l)/2 -1/2 
gxx grr 



u{r) {glx + B^) 



E. 



(3.10b) 



When 5 = 0, these reduce to the equations of motion in refs. [28, 33, 49], and when B is 
non-zero, these agree with the equations of motion in refs. [38, 50, 51]. Notice in particular 
that non-zero B leads to couplings between E and advertised. 

All of our results will follow from solutions of eqs. (3.10a) and (3.10b), with appropriate 
boundary conditions. Eqs. (3.10a) and (3.10b) are second-order, hence for each field, E 
and ay, we need two boundary conditions to specify a solution completely. At the future 
horizon, a solution for E or ay looks like a linear combination of in-going and out-going 
waves, with some normalizations. To obtain the retarded Green's function, we require that 
the fields be purely in-going [52-55]. For each field, the normalization then provides us 
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with a second boundary condition. Inserting in-going solutions into the action 3^2, we 
obtain G^{uj,k) via functional differentiation,^ 



G^^\uj,k) = lim (3.11) 
' ' r^oo 6a„6a* ^ ' 



Notice that 



G*i(a;, k) = lim = lim 

r— >oo oatOa^ r— s>oo 



5E 



5 at 



a? _ ,2ryEE 



6E6E* 



k'G%^{uj,k), (3.12) 



so any pole in G^^{u,k) not of the form 1/fc^ will produce a pole in G^^{uj,k). Given 
Gj^{uj,k), we obtain the spectral function x'^'^i'^^k) via its definition, eq. (3.2). 

Most of our results will be numerical. We will use the numerical method of ref. [57] 
to extract the locations of poles in G^(a;, A;) and the value of x^^i^i k), which we will now 
briefly review. Let us define a vector of fields 

Vir...k)^('"^^-'-'-'A. (3.13) 

yay{r,u,k) J 

We need four boundary conditions to specify a solution for V completely. As stated above, 
at the future horizon we choose solutions to be purely in-going waves. We next need 
to fix the normalizations of these in- going waves. Explicitly, the vector of near-horizon 
normalization factors, V^h, is, when T is non-zero,^ 

Vnh = ^hrn^ exp ^iw j dry/grr/\gtt\^ V{r,u}, k). (3.14) 

Notice that V^h is constant, independent of r, u, and k. On the right-hand-side of eq. (3.14), 
the exponential factor is designed to cancel the exponential factor that represents an in- 
going wave at the future horizon. 

In the method of ref. [57], we first pick any convenient value of V^h and solve the equa- 
tions, then we pick another V^h, linearly independent of the first, and solve the equations, 
thus constructing a basis of solutions. We can then write any solution in that basis. For 
example, suppose we solve the equations twice, with two normalizations V^^ = (1, 0)^ and 
^nh^ = (0,1)-^. Let the corresponding solutions be V^^^ and F^^^. Using these solutions, 
we define a 2 x 2 matrix P{r,uj, k) as 

P{r,u;,k) = (^y(i)(r,LJ,A;),y(2)(r^^^fc)^ ^ (3.I5) 



^Generically in AdS/CFT the on-shell action is divergent due to integration all the way to the AdS 
boundary. To compute renormalized correlation functions from the on-shell action, we first regulate the 
divergence by introducing a cutoff in r, introduce covariant counterterms localized on the cutoff surface, 
and then remove the regulator [56]. In our cases, the only near-boundary divergence that appears is 
oc F^vF'^" logr in the S^2 of the D7-brane, which we cancel with a counterterm. As discussed in ref. [10], 
a finite counterterm oc fp^F'"' is also present, with unfixed coefficient. As in ref. [10] we set the coefficient 
of the finite counterterm to zero. Our G'^{uj,k) in eq. (3.11) is then the renormalized Green's function. 

*We use non-zero T in all of our immerical analysis. In what follows, we consider T = only in 
section 5.1, where we compute the dispersion relation of holographic zero sound at T = 0. In that case the 
near-horizon behavior of V{r,Lj,k) is different from that at non-zero T: see eq. (5.4). 
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allowing us to write any solution as 

y(r,w,A;) =P(r,w,A;)Kh, (3.16) 

for some choice of V^h- Inserting a solution for V{r, uj, k) of the form in eq. (3.16) into 5^2 
and using eq. (3.11), we can express G^{oj, k) in terms of P{r,u}, k), 



G^'j^{oj,k) = -Mp lim r" 



(_^2 + ^2)-i 0^ 
1 



P'{r,uj,k)P~^{r,uj,k) 



(3.17) 



Generically, a pole will occur in G^{uj,k) when limr-s.00 det P~^(r, cj, A;) = 00, or 
equivalently when limr_s.oo det P{r, co, k) = 0. We can thus locate poles in Gj^{uj, k) simply 
by computing limj._s.oo det P{r, uj, k), rather than by computing the entire retarded Green's 
function, eq. (3.17). A pole in G^{oj,k) corresponds to a quasi-normal mode in the bulk, 
that is, a non-trivial solution for a gauge-invariant fluctuation with the in- going boundary 
condition at the horizon and a Dirichlet boundary condition at the AdS boundary [48]. 
Explicitly, imposing the Dirichlet condition at the AdS boundary means 

lim V{r, UJ, k) = lim P(r, uj, fc)Kh = 0. (3.18) 

r— >oo r— >oo 

If the solution V{r, uj, k), and hence V^h, is non-trivial, then eq. (3.18) is satisfied only when 
lim,.^oo det P{r, uj, k) = 0, demonstrating that a pole in G^{uj, k) is dual to a quasi-normal 
mode in the bulk. 

Let us summarize our numerical procedure. First, we choose values of T, d, B, uj, 
and k. Next we solve the equations of motion, eqs. (3.10a) and (3.10b), twice, with 
two linearly-independent values of V-ah- With the two resulting solutions we construct the 
matrix P{r,uj, k) as in eq. (3.15). At that point we can test for the presence of a pole in 
G^{uj,k) by testing whether lim^-^oo det P(r, tj. A:) = 0. We obtain the spectral function 
X^'^{uj, k) via eqs. (3.17) and (3.2). Repeating this process for various values of T, d, B, uj, 
and k, we obtained all of the numerical results that follow. 

4 Crossover at Zero Magnetic Field 

In this section we study the collisionless/hydrodynamic crossover in the D3/D5 system, in 
the limit B = The results are qualitatively similar to those in the D3/D7 theory [33]. 
We begin by reviewing the crossover in a LFL, which provides a useful reference point for 
studying the collective density excitations of our systems. 

4.1 Collective Modes in a Landau Fermi Liquid 

In a LFL at T = 0, the zero sound mode dispersion relation k{uj) at small uj/ (i has 

Re/c(a;)cxw, Im A;(a;) oc cj^. (4.1) 

If expressed instead as uj{k), then the dispersion relation at small k/ ^ has 

Rea;(A;) oc k, Ima;(A;) oc k^ . (4.2) 
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At non-zero T, the attenuation Im A; (a;) of the zero sound mode changes as T increases [2, 3]. 
The attenuation ImA;(a;) oc u, where is the quasiparticle colhsion frequency, 

1^ (X —- j^. (4.3) 

/i (1 + e-'^/^) 

On the right-hand-side of eq. (4.3), in the numerator the a;2 term is due to quantum interac- 
tions between the quasiparticles whereas the T'^ term is due to thermal scattering. At very 
low temperatures, T uj, the quantum term dominates and thus Imfc(ci;) oc uP' j \i. This 
is called the "coUisionless quantum" regime. As T increases, the thermal term begins to 
dominate when T ~ w, at which point the system enters the "coUisionless thermal" regime, 
where ImA;(w) oc T"^ j [i- At even higher temperatures, the collision frequency becomes sig- 
nificantly greater than the propagating frequency of the mode, u ^ oj. In this regime, the 
mode is very short-lived and so does not contribute significantly to the low-energy proper- 
ties of the system. However, in precisely this limit we expect hydrodynamic behaviour, such 
as hydrodynamic sound propagation and diffusion of charge, since the high frequency of 
thermal collisions brings about local thermal equilibrium. This collisionless/hydrodynamic 
crossover at ~ T'^/fJ' thus leads to a qualitative change in the response of the liquid to 
density perturbations. We summarize our discussion of the three regimes of a LFL with 
increasing T in fig. 2. 



(w/ju) 



1/2 



coUisionless quantum 
regime 



coUisionless 
thermal regime 



hydrodynamic 
regime 



Figure 2. The three regimes of charge density transport in a LFL, with increasing T//i. 



4.2 Collective Modes in Holographic Quantum Liquids 

We are interested in compressible states of the D3/Dp systems in the quantum liquid 
regime, that is, with a fixed, large density of matter present. We will thus normalise all 
dimensionful quantities by the appropriate power of the density, d, such that they are 
dimensionless. Specifically, we will define the dimensionless quantities 



T = 



irT 



UJ 



UJ 



k 



(4.4) 
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We will always work with T,iD,k <^ 1. Similarly, x^'^(a;, k) will denote the spectral function 
normalised to be dimensionless by an appropriate power of d. Notice that since T <C 1, 
normalising dimensionful quantities by powers of d is equivalent, up to numerical factors, 
to normalising by appropriate powers of the chemical potential /i, via eq. (2.10). 

Both of our holographic quantum liquids support a holographic zero sound excitation 
at T = [28]. For now, we will follow the convention of refs. [28, 29, 31-33] and write the 
dispersion relation of the holographic zero sound in the form uj{k), 



1 1 rl^/"- 

oj = ±^k-i±^P + 0(P) 
In 2n /io 



(4.5) 



(Recall that /io was defined in eq. (2.11).) These modes are called holographic zero sound 
because this T = dispersion relation has Reti;(/c) oc k and Imcj(/c) oc /c^, just like the zero 
sound excitation of a LFL, eq. (4.2) [28]. 

The holographic zero sound mode persists at small but non-zero temperatures in the 
range <^ T <^ k, with a dispersion relation described to a very good approximation by 
the T = form, eq. (4.5), as we show in fig. 3. This regime is analogous to the collisionless 
quantum regime of a LFL, where quantum interactions of quasiparticles dominate thermal 
collisions and hence the decay rate of zero sound — Ima;(A;) oc k'^. 



D3/D5: if =0,r =0.001 



D3/D5:B=0,r =0.001 




Figure 3. Left: Reel) as a function of k for the holographic zero sound dispersion relation in the 
D3/D5 system at S = and f = 0.001, within the range f^^f<t:k. Right: -Imw versus k 
for the same mode at the same T. In both figures the dots are our numerical results while the red 
lines come from the T — dispersion relation in eq. (4.5). Clearly the T — result in eq. (4.5) is a 
very good approximation to our numerical results for the low-temperature dispersion relation. 

As T increases and enters the range <^ k <^ T, the decay rate of the holographic 
zero sound mode increases as — Im£D(A:) oc T^, as shown in fig. 4. This is analogous to the 
collisionless thermal regime of a LFL, where thermal collisions of quasiparticles dominate 
quantum interactions, and the decay rate of zero sound —lmco{k) oc 

As T increases further still, the two holographic zero sound modes, viewed as poles of 
G^^{uj, k) in the complex Co plane, move down and also closer to the imaginary axis, following 
an approximately semi-circular path, and eventually collide on the imaginary axis to form 
two purely imaginary modes, as shown in fig. 5. The purely imaginary mode farther from 
the origin is short-lived, and does not contribute significantly to the low-energy physics of 
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Figure 4. Logarithm of the decay rate of the holographic zero sound mode in the D3/D5 system 
at -B = 0, normalized to the T = value, log (lmcj(T')/Im(I>(0)), as a function of logT for k — 0.01, 
within the regime ^ k. The dots are our numerical results, the vertical dashed black line denotes 
T — k, and the solid red best-fit line has gradient 2, indicating that lmu]{k) (x for T >k. 



the system. The purely imaginary mode closer to the origin is the hydrodynamic charge 
diffusion mode, with dispersion relation [49-51] 



3 11 

1' 2' 4' 



1 

2=14 



1. 



(4.6) 



as shown in fig. 5. The decay rate — Ima;(A;) of this hydrodynamic mode decreases as T 
increases. This is analogous to the hydrodynamic regime of a LFL. 
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Figure 5. Left: The complex Q plane, where the dots are our numerical results for the dominant 
poles of G*^{uj, k) in the D3/D5 system at _B = and k = 0.01 as T changes: the arrows point in 
the direction of increasing T. The two complex poles are due to holographic zero sound and the 
purely imaginary pole closest to the origin is due to charge diffusion. We define the location of 
the collisionless/hydrodynamic crossover as the value of T where the holographic zero sound poles 
collide on the imaginary axis. Right: — Imw(fc) versus k for the charge diffusion mode at T = 0.1, 
in the regime k <^T^ . The dots are our numerical results and the red line is the dispersion relation 
in eq. (4.6). 



These results for the D3/D5 system are qualitatively similar to those for the D3/D7 
system [33]. In both systems we define the location of the collisionless/hydrodynamic 
crossover as the value of T where the holographic zero sound poles collide on the imaginary 



-17- 



axis to form the purely imaginary poles. In both of the D3/Dp systems, we find that the 
crossover occurs when |a)| tx and k oc T^: for the D3/D5 system we show this in fig. 6. 
To find the exact proportionality constant in |w| cx T^, we performed numerical fits to the 



D3/D5:fl=0 D3/ 05:^=0 




-7 -6 -5 -4 -3 -2 _7 _6 _5 _4 _3 _2 

log ^ log r 



Figure 6. Left: log A: versus logT and Right: log|cZi| versus logT at the location of the coUi- 
sionless/hydrodynamic crossover for the D3/D5 system at _B = 0. hi both figures the dots are our 
numerical results and the solid red best-fit line has gradient 2, indicating that \oj\ cx and k oc 
at the crossover. 

data: we find that the crossover occurs when 

, , fo.30r2 forD3/D5, , , 

CD « < 4.7 
[0.45r2 forD3/D7. 

The crossover is also visible in the charge density spectral function xtt (Cd,k). In fig. 7 
we plot Xtt in the D3/D5 system at a fixed as a function of u), for various T. For 
small T, Xtt (t^) is dominated by a narrow peak due to the holographic zero sound mode. 
As T increases, the peak due to the holographic zero sound mode becomes wider as its 
decay rate —lmil}{k) increases: in fig. 5, the holographic zero sound poles move down 
into the complex uj plane. Near the collisionless/hydrodynamic crossover, the peak moves 
towards the origin which refiects the location of the corresponding poles of G^(a;, k) in the 
complex a) plane. As T increases further and the system enters the hydrodynamic regime, 
the peak near the origin becomes narrower due to the decreasing decay rate —lmoj(k) of 
the hydrodynamic charge diffusion mode. Again, these results for xtt ('^j^) in the D3/D5 
system are qualitatively similar to those of the D3/D7 system. 

In summary, in the D3/Dp systems, the three regimes of behaviour in response to 
density perturbations are very similar to those of a LFL. In particular, the decay rate 
—Im cjj(k) of the holographic zero sound mode has the same form as that of the LFL zero 
sound mode in both the collisionless quantum and collisionless thermal regimes. Addi- 
tionally, the boundaries between the regimes are at A; ~ T and A; ~ T^, so since |a)| ~ A; 
for the holographic zero sound mode (the poles in the complex u! plane in fig. 5 follow an 
approximately semi-circular path), the boundaries are located at cD ~ T and ui c^T^, as in 
a LFL. These similarities between the D3/Dp systems and LFLs are surprising because, as 
mentioned in sections 1 and 2, the D3/Dp systems are not LFLs. Indeed, no evidence has 



-18- 



D3/D5: B =0, * =0.01 D3/ D5: B =0, * =0.01 
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Figure 7. The charge density spectral function xtt in the D3/D5 system at S = as a function 
of ui and at fixed k = 0.01. Moving from the tallest peak to the shortest peak corresponds to raising 
the temperature: T = 0.01 (black), f = 0.032 (green), f = 0.1 (yellow), f = 0.16 (magenta), 
T — 0.32 (blue). The right-hand plot shows the tallest peak in full. 



been found of a Fermi surface in the D3/Dp systems, and in fact the low-energy dynamics 
may be governed by a (O-l-l)-dimensional CFT [11, 31, 32, 38]. 

In both of the D3/Dp systems, J^'s retarded Green's function, G^{u}, k), also exhibits 
poles in the complex O plane. As discussed in section 3, when B is non-zero and parity is 
broken, these poles will mix with those of G^^{uj, k), and hence will affect the crossover. As 
preparation for the non-zero- -B case, let us now discuss the poles in G^{uj, k) when B = 0. 
When B = 0, in both of the D3/Dp systems G^j^{u}, k) has no poles corresponding to long- 
lived modes, by which we mean poles satisfying a; — )• as A; — )• 0. The pole in G^^{u, k) 
closest to the origin of the complex ui plane is purely imaginary, but is not diffusive: for 
this mode ui (A;) is approximately independent of k. At low temperatures, for this mode 
lmu)(k) oc r^, as shown in fig. 8. 



D3/D5:i:=0.01 .B=0 D3/D7: A-=0.01, B=0 
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Figure 8. Left: log (— Im(Ii(fc)) as a function of log T for the pole in G^j^{uj, k) closest to the origin 
of the complex ui plane in the D3/D5 system at i? = 0, with k = 0.01, within the coUisionless 
regime k ^ T^. Right: The same as in the left figure, but in the D3/D7 system. In both figures 
the dots are our numerical results and the red best-fit lines have gradient 2, indicating that in both 
cases Imd)(fc) cx T^. 
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Following the convention of refs. [28, 29, 31-33], we described the poles in G*^{u),k) 
using complex frequency u>{k) with real momentum k. We could then define the crossover 
by a collision of poles in the complex Co plane, and could identify three distinct regimes of 
response to density perturbations, dependening upon the relative magnitudes of k, T and 
T^. For the holographic zero sound mode |a)| oc k, so we could equivalently characterize 
the three regimes by the relative magnitudes ofO, T and T^. Crucially, when B is non-zero 
such a simple translation from real k to real Co will in general not be possible, as we will 
see in the next section. 



5 Crossover at Non-zero Magnetic Field 

In this section we study the coUisionless/hydrodynamic crossover in both of the D3/Dp 
systems, with non-zero B. We will normalise B by the appropriate power of d to be 
dimensionless, 

We will always work with a large, fixed density, such that B,T,Cd,k <^ 1. In particular, 
we will always work with values of B below any known phase transitions. In general, the 
critical value of B for the known transitions is a complicated function of T, but at any given 
T the critical value of B is bounded from below by the T = values given in section 2. 
We will work with T > 10"'^ and values of B/T^ up to 10^ for the D3/D5 system and up 
to 10^ for the D3/D7 system, well below the known transitions. 

In section 5.1 we determine the dispersion relations of the holographic zero sound 
mode in our systems for T = and non-zero B analytically (without numerics), following 
refs. [38, 39]. In section 5.2 we determine the dispersion relation of the charge diffusion 
mode at non-zero B, following refs. [50, 51]. In section 5.3 we repeat the analysis of 
section 4 but with non-zero B, that is, we study the poles of G^^{uj, k) in the complex 
uj plane with non-zero B. For sufficiently small B we will show that a collision of poles 
occurs as we increase T, similar to the B = case, but that such a collision does not occur 
for sufficiently high B. We thus turn to the spectral function xtt ^) in section 5.4, 
which will show a clear crossover from holographic zero sound to charge diffusion as T 
increases with respect to u. The crossover is simplest to understand by studying the poles 
of G*^{uj, k) in the complex momentum plane, as we discuss in section 5.5. We summarize 
our results in section 5.6. 

5.1 Holographic Zero Sound at Non-zero Magnetic Field 

In this subsection we derive the dispersion relation of the holographic zero sound mode 
when T = and with B non-zero. For the D3/D5 system we reproduce the results of 
ref. [38], while for the D3/D7 system our results are novel. We will also show, using our 
numerical solutions, that the T = dispersion relation is a very good approximation to 
the low-temperature dispersion relation, as was the case when B = (see fig. 3). 

To determine the dispersion relation without numerics, we follow the procedure of 
refs. [28, 38]: we solve the equations of motion, eqs. (3.10a) and (3.10b), in two differ- 
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ent limits and then match the two solutions in a regime where the limits overlap. More 
specifically, we take a near-horizon limit of the equations of motion, meaning r — ?■ 0, solve 
the resulting equations, and then take a low-frequency limit of the solution. We then re- 
peat these operations in opposite order: we take a low-frequency limit of the equations of 
motion, solve the resulting equations, and then take a near-horizon limit of the solution. 

We begin by setting T = 0, or equivalently ru = 0, in eqs. (3.10a) and (3.10b), in which 
case \gtt\ = g^r — g^x = f"^- Our near- horizon r — ?• limit consists of taking r/d^/"^ <^ 1 
with r/y/B ^ 0(1), so that 

5^"' {gL + B') +<f^ d\ u? - uirfk^ ^ u?. (5.2) 

The resulting equations have solutions 

r / Ld u) J 

with constant normalization vector c, which is the T = analogue of the 14ih defined for 
non-zero T in eq. (3.14). On the right-hand-side of eq. (5.3) the e"^/'' factor describes an 
in-going wave as r — t- 0, as we discussed in section 3. 

We now want to take the low-frequency limit of the solution in eq. (5.3), which means 
w/r ^ 1 [28]. As emphasized in ref. [38], however, Kohn's theorem suggests that the 
holographic zero sound dispersion may exhibit a gap when B is non-zero, so if we wish to 
obtain a non-trivial result we must be careful not to take w small with fixed B. Instead, 
we should take a; small with B juj fixed, i.e. we should scale the gap with uj as we make u; 
small. In this limit the solution of eq. (5.3) becomes 



(r -|- iiS) 



y(r,a;,A;) = 7b . . 1 (5.4) 



We now need to perform the same operations in the opposite order. First, we take the 
low-frequency limit of eqs. (3.10a) and (3.10b), meaning cj/r ^ 1 and kjr keeping 
ujjk fixed. As above, when we take cj to be small we also take B to be small, in the sense 
that we keep Bju fixed. In these limits, in eqs. (3.10a) and (3.10b) we drop all terms not 
involving derivatives of E or a^, as these are sub- leading, and moreover in the remaining 
terms we drop all factors of B. The resulting equations have solutions 

fc2 -nu? (\ \ 13 1 
^— -2Fx 



E(r,u},k) = E^^^ + CEr'^~'^ , , . ^ 1 , , , , o 

^ ' ' ^ 1^1 + /r2" n(n- 1) ^2' 2 2n' 2 2n' r2« 

where E'^^\ ce, a'y \ and Cy are independent of r but can depend on uj and k. In an 



expansion near the AdS boundary, E^^"^ and a^^ are the leading, non-normalizable terms. 
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and hence act as sources for the dual operators and J^. Now we take the near- horizon 
hmit, r/d^" <C 1, 



E{r, CO, k) 
ay{r,uj,k) 



d 



UJ 



n 



+ CE-rr 
a 



,(0) 



d 



Cyd 



O 



O I r" 

,n+l ul+l/n" 



n+l 



'Id' 



(5.5a) 
(5.5b) 



We can now match solutions in the regimes where the two limits overlap. Specifically, 
we match the terms constant in r and the terms linear in r in eqs. (5.4) and (5.5). That 
gives us four equations involving the two components of c as well as E'^^\ ce, and Cy. 
We can use two of those equations to eliminate the two components of c, leaving us with 
two equations that we can express in matrix form: 



d ^ d 



UJ — — 



d 




(5.6) 



As discussed in section 3, for a quasi-normal mode the sources 

^(0) = and 4°^ = 0. We 

can then obtain a non-trivial solution for V{r, u, k) only if the matrix on the right-hand-side 
of eq. (5.6) has vanishing determinant. 

When B = 0, the matrix on the right-hand-side of eq. (5.6) becomes diagonal, i.e. 
E and ay decouple. In that case the determinant of the matrix will vanish when either 
diagonal entry vanishes. Setting the upper left entry to zero gives us the holographic zero 
sound dispersion relation at i? = 0, eq. (4.5). Setting the lower right entry to zero gives us 
a quasi-normal mode with dispersion u = —ifj-o, dual to a pole in G^{uj, k). This purely- 
imaginary, A;-independent mode is distinct from the purely-imaginary, fc-independent mode 
in G^^{uj, k) that we discussed near the end of section 4.2, which had Imcu(A;) oc (recall 
fig. 8), and hence ImcD(A!) — as T — )• 0. 

When B is non-zero, E and Uy couple, and so the poles in G^^{uj,k) and G^j^{uj,k) 
mix. Demanding that the determinant of the matrix on the right-hand-side of eq. (5.6) 
vanish when B is non-zero gives us 



UJ 



o- ^0 -3 



d^ 



n 



n 



0, 



(5.7) 



which determines the dispersion relation of the holographic zero sound: expressed as uj{k), 
when both B and k are non-zero and small, but of the same order, we find 



U!{k) 



-k^ ^ ^B^ - I 

n ^5 HQ 



1 

2n 



k' + 



B' 



+ O (P) . 



(5.8) 



When i? = we recover the result of ref. [28], namely eq. (4.5), and when B is non-zero, for 
n = 2 we recover the result of ref. [38]. In particular, we have found that the holographic 



zero sound dispersion relation is gapped: defining 



■'gap 



lim^_5^Q Reu;(A:), we have 



UJ. 



gap 



B. 



(5.9) 
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We can alternatively express the holographic zero sound dispersion relation as k(uj), 




(5.10) 



In fig. 9 we present our numerical results for the dispersion relation of the holographic 
zero sound in both D3/Dp systems with non-zero B, at low temperatures. We find that the 
T = dispersion relation, eq. (5.8), is a very good approximation to the low-temperature 
dispersion relation, as occured when B = (recall fig. 3). 



D3/ D5: r = 0.001 , B = 0.003 D3/ D5: T = 0.001 , fi = 0003 




D3/D7: r=0.001.fl=0.002 D3/D7: r=0.001, 5=0.002 




Figure 9. Top: Rea;(fc) (Left) and —lniu!{k) (Right) versus k for the holographic zero sound 
mode in the D3/D5 system with B ^ 0.003 and T = 0.001, in the range fc > T. Bottom: ReCj{k) 
(Left) and — Imw(fc) (Right) versus k for the holographic zero sound mode in the D3/D7 system 
with B = 0.002 and T = 0.001, in the range k ^ T. In all four plots the dots are our numerical 
results and the solid red lines come from the T = dispersion relation, eq. (5.8). Clearly the T = 
result eq. (5.8) is a very good approximation to the low-temperature dispersion relation, as was the 
case when B = (recall fig. 3). 



5.2 Charge Diffusion at Non-zero Magnetic Field 

As discussed in section 4.2, when B = a charge diffusion pole appears in G^^{uj,k) in 
the hydrodynamic regime k,Cd <^ . In this subsection, we will determine the dispersion 
relation of the longest-lived collective excitation in this hydrodynamic regime when B is 
non-zero. We will find that the longest-lived mode is that of charge diffusion. 
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Notice that previous holographic studies of hydrodynamics with a non-zero magnetic 
field, such as those of refs. [58-64], employed the typical "hydrodynamic" limit, in which 
T is larger than all other scales. In contrast, as mentioned at the beginning of this section, 
we work in a limit where the charge density is the largest scale. 

As in the previous subsection, to determine the dispersion relation of the charge dif- 
fusion mode without using numerics, we will solve the equations of motion, eqs. (3.10a) 
and (3.10b), in two different limits, the near-horizon and low-frequency limits, and then 
match the two solutions in a regime where the limits overlap. 

We begin with the near- horizon limit: we expand the equations of motion, eqs. (3.10a) 
and (3.10b), in (r — th) up to order (r — rn)- We then solve the resulting equations, and 
take the limits of low- frequency and low momentum, oj/T, k/T <^ 1, assuming oc \k'^\. 
In these limits we discard Oj^'s contribution to £"s equation of motion, being sub- leading. 
To find a quasi-normal mode, then, we only need the solution for E in these limits: 



E{r,uj,k) = I 1 + i 



ii.„ ("•^'■^)' 



9'uir)\ {9rr\r))\ 



ir-rn) ^nh, (5.11) 



where -Enh is the upper component of ^nh- Notice that here, unlike the previous subsection, 
we do not keep B/u fixed when we take the low- frequency limit. 

We now perform the same operations in the opposite order. First, in the equations 
of motion, eqs. (3.10a) and (3.10b), we take the low- frequency and low-momentum limits, 
w/T, k/T <^ 1, assuming \uj\ oc \k'^\. In those limits, in the equation of motion for E, 
eq. (3.10a), we discard the terms without derivatives on E or ay, as these are sub-leading, 
and moreover we take 

uj^ -uir fk"^ ^ -uirfk"^, (5.12) 

since we assume oc \k^\, which is suppressed relative to k"^ in our limits. The solution 
for E in these limits is then 



E{r, uj, k) 



E^^^+C df , _^ , ^^^^^ ' -j^, 5.13 



where E^^"^ and C are independent of r but can depend on cj and k. Next we perform the 
near-horizon limit, expanding the solution in eq. (5.13) in (r — ru) up to order (r — ru). 
We then match the constant term and the term linear in {r — ru) to the corresponding 
terms in eq. (5.11). The matching gives us two equation for E^^^ and C . Using one of those 
equations to eliminate C, and using the fact that each of guir), g~^{r), and u(r)^ has a 
simple zero ar r = rjj, we find 



" - (aUrn) + B^) V iuir^ry giV-'^/\r)\gu{r)\g;r'^'( 




(5.14) 
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Imposing the Dirichlet condition, E^'^^ = 0, we thus find the dispersion relation of a charge 
diffusion mode, 



uj{k) 



dr 



g'il^'^^^'^ (r ) I gtt {r) \gr}^'^ (r) 



In the D3/D5 system we can perform the integral in eq. (5.15), giving us 
_2 ^1 + 52 ^ / f2B^ + f^ + 2^2 (1 + f4) 



k^ + O {k^) . 

(5.15) 



Oik) 



+ 



T3 1 + 2^2 



;i + s2)' 



2-^1 




When S = we thus recover the result of eq. (4.6), and when B is non-zero we recover the 
result of refs. [50, 51]. In the D3/D7 system we have been unable to evaluate the integral 
in eq. (5.15) in closed form: 



-iT 



,Vl + B^T^ + T^ 
f^ + B"^ 



dp- 



;i + 52p2+p6): 



(5.17) 



When i? = we recover the result of refs. [65, 66], while to our knowledge our result with 
non-zero B is novel. 

While we have expressed the dispersion relations for the diffusive modes as uj{k), the 
expression for k{uj) is trivial to obtain by inverting eqs. (5.16) and (5.17). Notice that, 
when expressed as k{uj), the dispersion relation of the charge diffusion mode will have a 
non-zero real part. 

In fig. 10 we present our numerical results for the dispersion relation of the charge 
diffusion mode in both D3/Dp systems with non-zero B in the regime k <^ T^, which are 
extremely good approximations to the dispersion relations in eqs. (5.16) and (5.17), as was 
the case when B = (recall fig. 5). 

5.3 Poles in the Complex Frequency Plane 

Having found the dispersion relations of the long-lived excitations in both the T = limit 
and the hydrodynamic limit, we will now numerically investigate the collisionless/hydrodynamic 
crossover between these two limits. In this subsection we repeat the analysis of section 4.2, 
studying the motion of the poles in G^(a;, k) in the complex to plane as we increase T. 

Recall from section 4.2 that in both D3/Dp systems, when B = a collision of 
G^(a;,A;)'s poles occurs in the complex Cj plane: as we increase T at fixed k, the holo- 
graphic zero sound poles move (approximately) along semi-circles before colliding on the 
imaginary O axis to form two purely imaginary poles. The purely imaginary pole closest 
to the origin is the charge diffusion mode. Indeed, as T continues to increase, the charge 
diffusion pole approaches the origin, i.e. becomes more stable, while the other purely imag- 
inary pole moves away from the origin, becoming less stable. We summarize these B = 
results, for both D3/Dp systems, in fig. 11. 
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D3/D5: T = 0.32 , B = 0.003 D3/D7: 7=0.13, S=0.002 
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Figure 10. Left: — Imw versus k for the charge diffusion mode in the D3/D5 system with B = 
0.003 and f = 0.32, within the range k < T^. Right: The same for the D3/D7 system, with 
B = 0.002 and T — 0.13. For both modes. Re (w) = 0. In both plots the dots are our numerical 
results and the solid red lines are the dispersion relation in eq. (5.15). 




Figure 11. Left: The space (Retl), Imu), T), with blue dots indicating our numerical results for 
the dominant poles in 6*^(0;, A;) for the D3/D5 system with B — Q and k — 0.01. This plot is a 
three-dimensional representation of the left plot in fig. 5. Right: The same as the left plot, but for 
the D3/D7 system with B = and k = 0.01, and with black dots indicating our numerical results. 
In both cases we see that as we increase T (moving vertically in the plot) the two holographic zero 
sound poles move along (approximate) semi-circles and eventually collide on the imaginary w axis, 
producing two purely imaginary poles. The purely imaginary pole closest to the origin is that of 
charge diffusion, while the other purely imaginary pole is short-lived. 

When B is non-zero, the poles of G^^{ui, k) mix with those of k), and thus the 

latter has an "extra" purely imaginary pole near the origin, which we henceforth call the 
"transverse pole." For sufficiently small B a collision of poles still occurs, however. At 
very low T, the three poles of G^^ioj, k) closest to the origin are the two holographic zero 
sound poles and the transverse pole. As we increase T, the transverse pole moves down 
the Imo) axis, whereas the holographic zero sound poles again collide on the Ima) axis 
and produce two purely imaginary poles at a point closer to the origin than the transverse 
pole. One of these purely imaginary poles is the charge diffusion mode, which becomes more 
stable, moving up the Imw axis, as we increase T further. The second purely imaginary 



0.00003 
0.000025 
0.00002 
£ 0.000015 
0.00001 
5. X U)-^ 





0.0000 0.0002 0.0004 0.0006 0.(HH)8 0.0010 

I 



-26- 



mode becomes less stable as we increase T, and eventually collides with the transverse 
pole to form two propagating modes, i.e. two poles with nonzero real parts. Increasing 
T even more, these propagating modes move down into the complex O plane, becoming 
very short-lived, with decay rates — Ima)(A;) increasing as T increases. We present some 
small-i5 results in fig. 12, where we depict the motion of poles in the complex a) plane with 
B = 0.003 for the D3/D5 system and B = 0.002 for the D3/D7 system. 




Figure 12. Left: The space (Retlj, Imw, T), with blue dots indicating our numerical results for 
the dominant poles in G*^{ijj,k) for the D3/D5 system with B = 0.003 and k = 0.01. Right: The 
same as the left plot, but for the D3/D7 theory with B = 0.002 and k — 0.01, and with black 
dots indicating our numerical results. In both plots, as we increase T (moving vertically in the 
plot), the purely imaginary transverse pole moves down the hnui axis, while the holographic zero 
sound poles eventually collide on the Imw axis, above the transverse pole, producing two purely 
imaginary poles. One of these, corresponding to charge diffusion, moves up the Imui axis, while the 
other moves down, eventually intersecting the transverse pole and producing two propagating but 
short-lived modes (the "horseshoe" in the plot). 

We emphasise that at a fixed value of T, i.e. on a horizontal slice of figure 12, two kinds 
of pole appear near the origin of the complex uj plane. At low temperatures, these are the 
transverse pole and the poles of the holographic zero sound modes. At high temperatures, 
the three poles are the charge diffusion mode and the short-lived propagating modes. 

For small B we can define the location of the collisionless/hydrodynamic crossover by 
the collision of poles in the complex O plane, as in the B = case, however the story is 
not so simple for larger B. For sufficiently large B, as we increase T the holographic zero 
sound poles begin to move towards the Imw axis before turning around and remaining as 
propagating modes at higher T. Most importantly, the holographic zero sound poles never 
intersect the Ima) axis, i.e. no collision of poles occurs. As we increase T, the transverse 
pole moves down the ImO axis, reaches a minimum, and then moves back up the axis, 
eventually becoming the charge diffusion pole. We present some large-.B results in fig. 13, 
where we depict the motion of poles in the complex Co plane with B = 0.005 for the D3/D5 
system and B = 0.003 for the D3/D7 system. 

At even higher values of B than discussed above, the holographic zero sound poles 
do not appear to approach the ImtD axis at all, but rather seem to move directly away 



-27- 



Figure 13. Left: The space (ReLu,lmu!,T), with blue dots indicating our numerical results for 
the dominant poles in G'^^{uj,k) for the D3/D5 system with B = 0.005 and k = 0.01. Right: The 
same as the left plot, but for the D3/D7 theory with B = 0.003 and k — 0.01, and with black dots 
indicating our numerical results. In both plots, as we increase T (moving vertically in the plot), 
the holographic zero sound poles move towards the Ima) axis, reach a minimum value of Red), and 
then move away from the ImcD axis, remaining as propagating modes at higher values of T. The 
purely imaginary transverse pole moves down the Imw axis, reaching a minimum value, and then 
moves back up the axis, eventually becoming the charge diffusion pole. 

from it as T increases. The transverse pole's behavior is similar to that in fig. 13, i.e. as 
T increases the transverse pole moves down the Ima) axis, reaches a minimum, and then 
moves back up the axis, eventually becoming the charge diffusion mode.^ 

As at smaller B, at a fixed temperature T in the large- -B regime, i.e. on a horizon- 
tal slice of fig. 13, two kinds of mode are always present, one purely imaginary and the 
other propagating, however, unlike the small-S regime, here we can show that the purely 
imaginary mode is charge diffusion and the propagating modes are holographic zero sound. 
Consider for example the temperature range k ^ T where, when B = 0, the holographic 
zero sound mode is present, with a dispersion relation that approximately obeys the T = 
dispersion relation eq.(4.5), as shown in fig. 3. In fig. 14 we present our numerical results 
for the dispersion relation of the holographic zero sound as a function of B in both D3/Dp 
systems. We find that our numerical results are very well approximated by the 7^ = 
result, eq. (5.8), for all B. In fig. 15, we present our numerical results for the dispersion 
relation of the longest-lived purely imaginary pole as a function of B in the temperature 
range k ^ T. We find that for large enough B, this purely imaginary pole is simply the 
charge diffusion mode. In short, we find that for a fixed S> T, both the holographic zero 
sound and charge diffusion modes coexist at sufficiently high B. 

^We have presented numerical results for poles of G%{lj, k) in the complex lo plane with fixed momentum 
k and increasing T. If we instead fix T and increase k, then we find results similar to those of the 
compressible states studied, using holography, in ref. [41], namely A/" = 4 SYM coupled to (2+l)-dimensional 
fiavor fermions alone, rather than a hypermultiplet, with non-zero U{l)b charge density. In that case, the 
holographic dual was a probe D-brane with non-zero worldvolume electric and magnetic fiuxes, including 
magnetic fluxes along directions of the internal space wrapped by the D-brane. 
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Figure 14. Top: Re a; (Left) and — Imw (Right) as functions of B for the holographic zero sound 
poles of the D3/D5 system with k = 0.01 and f = 0.001, within the range A: » T. Bottom: The 
same as the top plots, but for the D3/D7 system. In all plots the dots are our numerical results and 
the sohd red lines come from the T = dispersion relation in eq. (5.8). Clearly eq. (5.8) provides 
a very good approximation to the numerical results. 



D3/D5:r=a001 ,«:=0.01 



D3/D7: r=0.Q01,t-=0.01 
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Figure 15. Left: — Imw as a function of B for the longest-lived purely imaginary pole of the 
D3/D5 system with k — 0.01 and T = 0.001, within the range k ^ T. Right: The same as the 
left plot, but for the D3/D7 system. In each plot the dots are our numerical results and the solid 
red line comes from the dispersion relations in eq. (5.16) or (5.17). For sufficiently large B, we see 
that the purely imaginary pole becomes that of charge diffusion. 



In physical terms, we have learned that if we perturb these systems at a fixed mo- 
mentum k, then for sufficiently large B we will see both holographic zero sound and 
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charge diffusion, hence for large B we cannot characterize the reponse of the system as 
either colhsionless or hydrodynamic. We then seem to have no way of defining a cohision- 
less/hydrodynamic crossover. The story wih be very different if we perturb these systems 
with fixed real frequency cu, as we will now show. 

5.4 Charge Density Spectral Functions 

The poles of G*^{uj, k) in the complex O plane are not directly observable. To gain per- 
spective, we thus turn to a quantity that is directly observable (in principle), the charge 
density spectral function xtt (w, fe). 

We have so far considered perturbations with fixed real momentum k. In such cases, 
the charge density spectral function xtt as a function of Co, will exhibit peaks at 

frequencies u determined by the positions (and residues) of the poles in the complex cu 
plane. For example, when B = 0, as we increase T, the spectral function always exhibits 
only a single dominant peak that moves towards the origin, reflecting the movement of the 
corresponding pole in the complex ui plane towards the Imw axis, as we saw in fig. 7. 

In fig. 16 we present our numerical results for logxtt ('^, A;) in the (logc<j/T^,logA;/T'^) 
plane for both D3/Dj) systems at two different values of B/T^. Here we chose to normalize 
(D and k by because in the B = case Co/T'^ and k/T"^ are order one at the crossover. 
Given those normalizations and the fact that the gap in the T = holographic zero sound 
dispersion relation goes like Wgap oc B, the natural normalization for B was then B/T"^. 
To make the relationship between poles of G^^{LJ,k) and peaks in xtt (w,A;) explicit, in 
fig. 16 we also superimposed on the plots our numerical results for the locations of poles 
of G^^{uj, k) in the complex a) plane (the red dots and blue crosses) as well as the absolute 
values of the dispersion relations of the long-lived modes, i.e. we plot the value of |a) (A!)| 
for the holographic zero sound dispersion relation in eq. (5.8) and the charge diffusion 
dispersion relations in eqs. (5.16) and (5.17).^ 

Consider first the plots of logXtt (i^, fc) with B/T"^ <^ 1, figs. 16 (a) and (c). For a 
fixed value of the momentum k, i.e. along a horizontal line across the contour plot, the 
spectral function exhibits only a single peak, which for k ^ T^, near the top of each plot, 
is due to the holographic zero sound pole, as indicated by the blue crosses and solid blue 
line, and for k <^ T^, near the bottom of each plot, is due to the charge diffusion mode, as 
indicated by the red dots and solid red line. When A; ~ T^, these two peaks join, refiecting 
the collision of the corresponding poles of G^^{uj,k) in the complex ui plane. The "extra" 
poles due to mixing with G^^{oj,k) do not produce any significant peaks in the spectral 
function at these low values of B/T"^. These plots provide additional evidence that at low 
B and with fixed k, a crossover occurs similar to the crossover at = 0. 

Now consider the plots of logxtt (w, A;) with B/T'^ ^> 1, figs. 16 (b) and (d). These are 
qualitatively different from the B/T^ <C 1 plots: now, at a fixed value of k (a horizontal 
line across fig. 16 (b) or (d)), the spectral function exhibits two peaks, not one peak. At 

®In general, the precise location of the peak in xtt (<^, k) actually depends on the location and residue of 
the pole in Gj|(cj, k) in a complicated way. The results depicted in fig. 16 show that for our systems \Gj (fc) | 
gives a very good approximation to the location of the peak. Notice that for a long-lived, propagating 
mode, this agrees with the usual expectation that the peak is located at the propagating frequency. 
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(a) D3/D5: logxtt at T = 10"^/^ B/T^ = 10"^ (b) D3/D5: logxtt at T = IQ-^, B/T^ = 10^ 
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(c) D3/D7: logoff at T = 10-\B/f^ = 10"^ (d) D3/D7: logxtt at T = 10'^, B/f"^ = 10'' 
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Figure 16. (a) Contour plot of our numerical results for logxtt(w, k) over the (loga;/T2,logfc/T2) 
plane in the D3/D5 system with T = 10^^/^ and B/T^ = 10~^. Lighter colours represent larger 
values of the spectral function, with logXft — — 15, — 12.5, — 10, . . . , 2.5 from darkest to lightest. 
The red dots and blue crosses are our numerical results for the locations of purely imaginary and 
propagating poles in G^(w, fc), respectively. The solid red and blue lines are the values of \ui{k)\ 
for the dispersion relations of the charge diffusion, eq. (5.15), and holographic zero sound, eq. (5.8), 
respectively, (b) The same as (a), but with T = 10"'^ and B/T'^ = . (c) The same as (a) and (b), 
but for the D3/D7 system with f = IQ'^ and B/f^ = 10"^ and with log Xtt = -10, -8, -6, . . . , 2 
from darkest to lightest, (d) The same as (c), but with f = IQ-^ and B/f"^ = 10''. 



very high k/T'^, near the top of each plot, the holographic zero sound peak is dominant, 
and at very low Ai/T^, near the bottom of each plot, the charge diffusion peak is dominant. 
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In the intermediate region, however, the holographic zero sound and charge diffusion peaks 
coexist with roughly equal spectral weight. These results agree with what we found by 
studying poles in the complex frequency plane in the last subsection. In particular, in 
figs. 16 (b) and (d) the holographic zero sound and charge diffusion peaks are no longer 
joined: as we saw in the previous subsection, the corresponding poles in the u) plane do 
not collide. 

We have so far merely confirmed the analysis of the previous subsection. Now, however, 
let us consider Xttioj, k) as a function of A; at a fixed frequency Co, i.e. a vertical line in any 
of the plots of fig. 16. With fixed a), we find that for any B and any a) the spectral 
function exhibits only one peak, which at low w/T^, near the left in each plot, is due to 
charge diffusion, and at high lo/T'^, near the right in each plot, is due to holographic zero 
sound. We also find that for larger B, figs. 16 (b) and (d), at intermediate frequencies 
the spectral function is suppressed by orders of magnitude relative to the holographic zero 
sound and charge diffusion peaks. That suppression is due to the gap in the holographic 
zero sound dispersion, as we discuss in more detail in the next subsection. At smaller 
B, such as figs. 16 (a) and (c), our numerical results indicate that if some suppression of 
spectral weight occurs in some non-zero range of frequencies at small B, then either the 
suppression is very small or the range of frequencies is very narrow. 

By studying the spectral function xttii^,k), we have learned that the response of the 
system is simplest to understand using perturbations with fixed frequency Co rather than 
fixed k: with fixed Co the spectral function always exhibits only a single peak for any 
B, in which case we should be able to characterize the response of the system as either 
collisionless or hydrodynamic. We are thus able to define a collisionless/hydrodynamic 
crossover, as we show in the next subsection. 

5.5 Poles in the Complex Momentum Plane 

We saw in the previous subsection that the spectral function xtti^^k), as a function of 
momentum A; at a fixed frequency Co, exhibits only a single peak, which in physical terms 
means the response of the system to an external perturbation with frequency Co is dominated 
by a single excitation. To study that excitation, we will determine the poles of G^^{uo, k) in 
the complex k plane. We will look for poles with Re(A;) > and lva.{k) > 0, and will focus 
on the poles closest to the origin of the complex k plane, corresponding to the modes with 
the smallest propagating momentum Re(A;) and attenuation Im(A;). 

First, let us revisit the known collisionless/hydrodynamic crossover at S = 0, now 
studying the behavior of poles in the complex k plane. In fig. 17 we plot log (Re k {Co) /T^) 
and log (imfe {Co) /T^) of the pole closest to the origin as functions of log (a)/T^), for both 
D3/Dp systems with B = Q and T = 0.001. In fig. 17 we also plot the dispersion relations 
of the holographic zero sound mode at T = 0, eq (5.10), and of the charge diffusion poles 
in the complex k plane, obtained by inverting eqs. (5.16) and (5.17), as well as the location 
of the collisionless/hydrodynamic crossover extracted from the collision of poles in the 
complex Co plane, eq. (4.7). We see that at the crossover, no collision of poles occurs in the 
complex k plane. Instead, we see a single pole whose dispersion relation k{Co) changes at 
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the crossover. For example, we can see from the left plots of fig. 17 that ReA;(a;) exhibits 
a "kink" at the crossover in both D3/Dp systems. 



D3/D5;B =0,r =0.001 D3/D5:B =0,r =0.001 




-5 5 10 -5 5 10 

Log(iJ/f') LogWf') 



Figure 17. Top: log (Re fc (tl>) /T^) (Left) and log (im fc (w) /T^) (Right) as functions of 
log(w/T'^) for the pole in G*^{uj,k) closest to the origin of the complex k plane in the D3/D5 
system with B — and T = 0.001. The dots are our numerical results, the solid red lines, be- 
ginning on the left of each plot, come from the dispersion relation of the charge diffusion mode 
obtained by inverting eq. (5.15), the solid blue lines, beginning on the right of each plot, come from 
the holographic zero sound dispersion relation at T = 0, eq. (5.10), and the vertical dashed line 
is the location of the coUisionless/hydrodynamic crossover, as defined by a collision of poles in the 
complex Co plane, eq. (4.7). Bottom: The same as the top plots, but for the D3/D7 system. On 
the right plots we also indicate the hydrodynamic, collisionless thermal, and collisionless quantum 
regimes, as defined in section 4.2. 

When B is non-zero, we again find no collision of poles in the complex k plane. In 
fact, when B <^T^, the behavior of the poles is very similar to the ^ = case, as we show 
in fig. 18, which is the same as fig. 17, but now with B/T'^ = 10~^. In particular, we still 
see a kink in Rek{uj) at approximately the same value of a)/T^ as at .B = 0. 

When B is non-zero and large, B ^ T^, we find a qualitative difference from the 
B <^ case: at low lj/T'^ the mode is that of charge diffusion until a local maximum 
in Iiek{id) at some ui. For larger O, the mode is no longer that of charge diffusion. We 
show this in fig. 19, which is the same as figs. 17 and 18, but now with B/T'^ = 10^ for the 
D3/D5 system and B/T"^ = 10^ for the D3/D7 system. Remarkably, the local maximum in 
Refc(a;) occurs at approximately the same value of oj/T"^ as the kink in ReA;(ti;) at .B = 
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D3/D5:r =0.1,B/r^ =10"- D3/D5: T =0.1 ,B / =10"^ 
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Figure 18. Top: log (Re fc (tZi) /T^) (Left) and log (im fc (w) /T^) (Right) as functions of 
log(a}/T'^) for the pole in 6*^(0;, fc) closest to the origin of the complex k plane in the D3/D5 
system with T = 0.1 and B/T^ — 10^^, within the regime B/T^ ^ 1. The dots are our numerical 
results, the solid red lines, beginning on the left of each plot, come from the dispersion relation of 
the charge diffusion mode obtained by inverting eq. (5.15), the solid blue lines, beginning on the 
right of each plot, come from the holographic zero sound dispersion relation at T = 0, eq. (5.10), 
and the vertical dashed line is the location of the coUisionless/hydrodynamic crossover, as defined 
by a collision of poles in the complex Q plane at B = 0, eq. (4.7). Bottom: The same as the top 
plots, but for the D3/D7 system. All four plots, in this small B regime, are qualitatively similar to 
the corresponding plots at B = in fig. 17. 

(the vertical dashed line in the left plots of fig. 19). 

In the intermediate regime B ~ T^, at low Cj/T"^ the mode is still that of charge 
diffusion while at high Cj/T'^ the mode is still that of holographic zero sound, and the 
transition between the two still occurs at approximately the same value of lj/T^ as the 
kink in Re k{u)) at i? = 0, however at the transition neither a sharp kink nor a clear 
maximum in ReA;(ci)) were visible to us. The spectral function, at fixed w, as a function of 
k, still exhibits only a single peak, which is due to charge diffusion for small u)/T'^ and to 
holographic zero sound for large w/T^. 

We have found that as we increase a) the pole in G^(a;, k) in the complex k plane 
ceases to be that of charge diffusion approximately at the location of the crossover in the 
B = Q case, eq. (4.7), for any value of B. In other words, for any value of S, the response 
of the system to a density perturbation with fixed a) below this value is dominated by 
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D3/D5: T = 10"', fi / r" = 10' D3/D5: T = 10-^,8/ = lO' 
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Figure 19. Top: log (Re fc (cj) /T^) (Left) and log (im fc (w) /T^) (Right) as functions of 
log(w/T^) for the pole in G^(a;,fc) closest to the origin of the complex k plane in the D3/D5 
system with T — 10^^ and B/T^ — 10'^, within the regime B/T^ ^ 1. The dots are our numerical 
results, the solid red lines, beginning on the left of each plot, come from the dispersion relation of 
the charge diffusion mode obtained by inverting eq. (5.15), the solid blue lines, beginning on the 
right of each plot, come from the holographic zero sound dispersion relation at T = 0, eq. (5.10), 
and the vertical dashed line is the location of the coUisionless/hydrodynamic crossover, as defined 
by a collision of poles in the complex oj plane at i? = 0, eq. (4.7). Bottom: The same as the top 
plots, but for the D3/D7 system with B/f'^ = 10^. 

charge diffusion, while the response to a perturbation with larger ui is not. We have thus 
identified a clean boundary to the regime of hydrodynamic reponse that is approximately 
S-independent. We therefore define the location of the crossover from the hydrodynamic to 
the collisionless regime as occuring at this value of w/T^, eq. (4.7): for the D3/D5 system 
ud ^ O.SOT^ and for the D3/D7 system Co ^ OAbf^. 

Let us now discuss in detail the region of suppressed spectral weight that we mentioned 
at the end of section 5.4, in the regime of large magnetic field, B ^ T^. In fig. 20 we 
reproduce fig. 16, the contour plots of logxtt(a;,fe) in the plane of (logw/T^, log fe/T^) , 
but now with two differences: first, instead of a contour plot of log xtt('^, A;) we present a 
density plot of xtti^,k), and second we superimpose on the density plots our numerical 
results for the locations of the pole of Gj^(a;, k) closest to the origin in the complex k plane 
as well as the values of |A:(a;)| for the T = holographic zero sound dispersion relation, 
eq. (5.10), and for the charge diffusion, obtained by inverting eq. (5.15). In fig. 20 we see 
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clearly that at a fixed frequency u the peak in the spectral function xtt (w, k) is given to 
good approximation by these values of \k {oj)\. 



(a) D3/D5: Xtt at T = 10"^/^ B/T^ = 10"^ (b) D3/D5: xtt at T = 10-^ B/T^ = 10^ 




Figure 20. (a) Density plot of our numerical results for Xtt{^,k) over the (logw/T^, log fc/T^) 
plane in the D3/D5 system with T = 10~^/^ and B/T^ — IQ^^. (This plot uses the same data as 
that used to make fig. 16.) Lighter colours represent larger values of the spectral function, as in 
fig. 16. The black dots are our numerical results for the location of the pole in G^^{u), k) closest 
to the origin in the complex k plane. The solid red line starting at the bottom left of the plot is 
the value of | k (ui) | for the charge diffusion mode obtained by inverting the dispersion relation in 
eq. (5.15), and the solid blue line starting at the top right of the plot is the value of \k {uj)\ for the 
holographic zero sound mode, eq. (5.10). (b) The same as (a), but with T = lO^'^ and B/T"^ ~ IQ^ . 
(c) The same as (a) and (b), but for the D3/D7 system with f = 10"^ and B/T^ = 10"^^ (d) 
The same as (c), but with T = 10"^ and B/T^ = 10^. 
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In the regime B/T"^ S> 1, as shown in figs. 20 (b) and (d), suppose we start with small 
/T^, that is, we consider a vertical line near the left of each plot, where the spectral 
function exhibits a single peak due to charge diffusion. As we increase w/T^, moving to 
the right in each plot, we encounter the collisionless/hydrodynamic crossover, where the 
peak is no longer due to charge diffusion. For larger a)/T'^, a peak is still present in the 
spectral function, but is smaller by orders of magnitude compared to the charge diffusion 
peak. Only at sufficiently large w/T^ does the peak grow again by orders of magnitude, 
now being the peak due to holographic zero sound. In other words, for an intermediate 
range of frequencies (D/T^, the spectral weight is suppressed. The high-frequency end of 
that intermediate region (on the right in the plots), is marked by a minimum in the value 
of \k ((D) I of the corresponding pole: in each of figs. 20 (b) and (d), the solid blue line on 
the right has a cusp with a distinct minimum. The value of u) at that minimum, (Dmim is 
fixed by the gap in the holographic zero sound dispersion relation, cDgap in eq. (5.9): in 
fig. 21 we plot log (B/r^) versus the value of log (wmm/T"^) j and find very good agreement 
with the value of Wgap. Notice that such a minimum does not occur for B/r2 < 1, i.e. the 
solid blue line in figs. 20 (a) and (c) has no minimum like that in figs. 20 (b) and (d). 




4 6 8 10 2 4 6 8 10 



Figure 21. Left: log (i?/T^) versus log (wmin/T'^) for the D3/D5 system, where Wmin is the high- 
frequency end of the region of suppressed spectral weight, i.e. the cusp in the solid blue lines in 
figs. 20 (b) and (d). The dots are our numerical results and the solid red line is derived from the 
gap in the T = holographic zero sound dispersion relation, the Wgap in eq. (5.9). Clearly eq. (5.9) 
is a very good approximation to the numerical result for Wmin- Right: The same as the left figure, 
but for the D3/D7 system. 

5.6 Overview: The Collisionless/Hydrodynamic Crossover 

We can now give an overview of the collisionless/hydrodynamic crossover in our systems 
at non-zero B. When B <^ T^, if we fix a real k and study poles in G^^{co,k) in the 
complex id plane, then we find that as we increase T a collision of poles occurs, allowing 
us to define the crossover in a fashion similar to the ^ = case. When B ^ T^, no such 
collision of poles occurs in the complex w plane. Indeed, in this regime of large B, the 
spectral function, for fixed A; as a function of w, generically exhibits two peaks, one for 
holographic zero sound and one for charge diffusion, preventing us from identifying this 
regime as either collisionless or hydrodynamic. We found that for any B, if we instead 
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fix real w and study poles in G^^{uj, k) in the complex k plane, then we find a single pole, 
and in the spectral function we find a single peak. As we increase T, we could then define 
the crossover as the value of T where the single pole becomes that of charge diffusion, and 
hence the system enters the hydrodynamic regime. Remarkably, we found that the value of 
T where that crossover occurs is approximately B-independent, being given by the 5 = 
value, eq. (4.7): ^ 0.30 for the D3/D5 system and |cj| 0.45 for the D3/D7 
system. In the language of LFL theory, this would indicate that the collision frequency 
of the quasiparticles is independent of B. We summarize our results in fig. 22, which we 
represented schematically in fig. 1. 




Figure 22. An overview of our main results (see also fig. 1): log (^B/T^) versus log (^uj/T^) for the 
D3/D5 system (Left) and the D3/D7 system (Right). The lines and dots have the same meanings 
in the two plots. The solid red vertical line denotes the coUisionless/hydrodynamic crossover, defined 
for increasing T with fixed ui as the value of T where the pole of G^ (a;, k) closest of the origin of the 
complex k plane becomes that of charge diffusion. That value of T is approximately the same as the 
B = value in eq. (4.7). The dots on the solid red vertical line above the dashed lines are extracted 
from poles in the complex k plane, as discussed in section 5.5, while those below the dashed lines 
are extracted from a collision of poles in the complex uj plane, as discussed in sections 4.2 and 5.3. 
In the region between the horizontal dashed lines we had difficulty extracting a precise location of 
the crossover from our numerics. The solid red diagonal line is the value of the gap in the T = 
dispersion relation for holographic zero sound, eq. (5.9). The dots on the solid red diagonal line are 
our numerical results for Wmin, reproduced from fig. 21. On the right figure we have indicated the 
nature of the dominant peak in the spectral function with fixed w as a function of k. 



We can actually provide a quantitative estimate of the 0(1) value of B/T'^ at which, 
as B increases in the collisionless regime, the gap in the holographic zero sound dispersion 
relation produces suppression of spectral weight. In fig. 22, the diagonal red line represents 
the value of the gap Wgap in the T = dispersion relation of the holographic zero sound, 
eq. (5.9), while the vertical red line represents the location of the crossover, i.e. the value 
of ui in eq. (4.7). By extending the diagonal red line down to the intersection with the 
vertical red line, we can estimate the 0(1) value of B/T"^ that we want: we simply equate 
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the Wgap oc B in eq. (5.9) with the lo ccT'^ 



in eq. (4.7) and solve for B/T'^, 



0.56 for D3/D5, 
0.63 for D3/D7. 



(5.18) 



6 Summary and Outlook 

We used hofography to study compressible states of the D3/D5 and D3/D7 theories that 
were not solids, superfluids, LFLs, or NFLs. In particular, we studied the effect of a 
non-zero magnetic field B on the collisionless/hydrodynamic crossover in these systems 
as we increased the temperature T, for values of B below any known phase transitions. 
We found that the crossover was simplest to understand by studying the spectrum of 
collective excitations of these systems with fixed frequency uj, where we could define the 
crossover for any value of B as the value of T where the dominant pole in G^^{uj, k) becomes 
that of charge diffusion. Remarkably, that value of T is, to a very good approximation, 
independent of -B. In the high-frequency, coUisionless regime, for values of B above a value 
B/T'^ ~ 0{1) estimated in eq. (5.18), we also saw a suppression of spectral weight due to 
the magnetically-induced gap in the holographic zero sound dispersion relation. Our main 
results are summarized in figs. 1 and 22. We believe that, given our minimal ingredients, 
our fig. 1 may be generic to compressible states described holographically by probe DBI 
actions with non-zero electric and magnetic flux. In other words, we believe that fig. 1 may 
be characteristic of such compressible states, and hence may help in classifying them. 

Various open questions remain about these, and other similar, holographic quantum 
liquids. Here we will list only a few such questions. 

In the D3/D7 system we only studied excitations with momenta perpendicular to 
the magnetic field, but we have another option: the momentum can also have a non- 
zero component parallel to the magnetic field. A natural question is how the non-zero 
magnetic field affects modes with momentum parallel to the magnetic field, for example, 
how does non-zero B affect the holographic zero sound mode that propagates parallel to 
the magnetic field? What happens to modes parallel to the magnetic field during the 
coUisionless /hydrodynamic crossover? 

We only considered values of B below any known phase transitions. What happens at 
higher values of B, in the symmetry-broken phases? Clearly the spectrum of excitations 
will then include the Goldstone bosons associated with the symmetry breaking. How do 
those affect the collisionless/hydrodynamic crossover? 

More generally, an open question is whether these compressible states of the D3/Dp 
theories are the true ground states, that is, whether these states are global minima of the 
free energy. Indeed, as we mentioned in section 2, for the D3/D7 system, at one point 
in the phase diagram in the grand canonical ensemble with non-zero B, these states are 

^For a holographic quantum hquid in which the U {l)b density is produced by fermions alone, a result of 
ref. [41] is that the holographic zero sound's dispersion relation becomes gapped when B/T^ > 0.18. The 
relationship between that value of B/T^ and the value above which the spectral weight of the holographic 
zero sound is suppressed would be interesting to study. 
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known to exhibit an instability towards an inhomogenous state. We hasten to add that 
even if the states we studied are not the true ground states, they are a useful theoretical 
laboratory that may reveal some guiding principle(s) to classify compressible states. 

Another perspective on the collective excitations of these systems in the collisionless 
regime, which we have not explored here, is given by the semi-holographic description [11, 
31, 32, 38], an effective description in which some critical degrees of freedom, described by 
a (0 + l)-dimensional CFT dual holographically to degrees of freedom in a near-horizon 
AdS2, produce the dissipation of the holographic zero sound mode. Clearly in our systems 
the bulk geometry has no near- horizon AdS2 (unlike extremal Reissner-Nordstrom-AdS), 
nevertheless quite generically fluctuations of the Dp-brane worldvolume fields effectively 
"see" AdS2 deep in the bulk, due to the worldvolume electric flux [11, 31, 32, 38]. An 
open question is what aspects of the collisionless/hydrodynamic crossover, at either zero 
or non-zero B, can be understood using semi- holographic techniques. 

The behavior of collective excitations with increasing T for the compressible states dual 
to Reissner-Nordstrom-AdS is qualitatively different from that of the D3/Dp systems [67, 
68]. Although those states do exhibit a sound mode in the collisionless regime [67], as T 
increases no collision of poles in the complex uj plane occurs [68]. Does a non-zero magnetic 
field produce a gap in the sound dispersion relation in those states? If so, how does that 
affect the behavior of the collective modes as T increases? 
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